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Abstract 

We investigate quantum gravity in the path integral formulation using the Regge 
calculus. Restricting the quadratic link lengths of the originally triangular lattice 
the path integral can be transformed to the partition function of a spin system with 
higher couplings on a Kagome lattice. Various measures acting as external field were 
considered. Extensions to matter fields and higher dimensions are discussed. 

More than thirty years ago Regge proposed to approximate Riemannian manifolds by piece- 
wise flat, simplicial manifolds. He developed a discrete description of General Relativity 
in which space-time is triangulated by a simplicial lattice, the Regge skeleton [1]. Thus, 
the lattice becomes a dynamical object, with the edge lengths describing the evolution of 
space-time. More precisely, dynamics comes from the variation of the link lengths keeping 
the connectivity of the triangulation fixed. We briefly mention that a somehow comple- 
mentary method is provided by the dynamical triangulation approach in which the length 
of the links is kept fixed and the incidence matrix is varied. 

Within the Regge scheme the Einstein-Hilbert action for lattice gravity in two dimen- 
sions is given by [2] 

S = \J2A t -2kJ2$v , (1) 

t v 

with A a cosmological constant and A t the area of triangle t. The curvature of the surface 
is concentrated at the vertices of the lattice. The deficit angle at a single vertex v is defined 
by 

8 V = 2tt — [vertex angle at v] , (2) 

tDv 

where tDv denotes the triangles meeting at v. 
The corresponding continuum action is 

S c = j cPxy/g(\ - kR) , (3) 

in which g is the determinant of the metric tensor and R the curvature scalar. According to 
the Gauss-Bonnet theorem the Einstein part of the action is a topological invariant equal 
to 47r times the Euler characteristic of the surface. Therefore, for a manifold with fixed 
topology the term proportional to k can be dropped. In the following only lattices with the 
topology of a torus will be considered. 
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In the path integral formulation a quantization of the above action proceeds by evalu- 
ating the expression 

Z = J dn[q]e- S[q] . (4) 

Unfortunately, a unique prescription for the measure does not exist, however, an appropriate 
choice proposed in the literature is [3] 



with m G JR defining a one-parameter family. 

The central idea of this investigation is to transform the path integral to a partition 
function of a spin system. Therefore, we allow all the squared edge lengths to take two 
values 

qi = 1 + tax , < e < 0.6, <r,eZ 2 . (6) 

The real parameter e is restricted to fulfill the Euclidean triangle inequalities for the q^s. 
To rewrite the action in terms of o\ we consider a single triangle (cf. Figure 1). 




Fig.l: Notation for our triangular lattice. Triangle / consists of edges with the quadratic lengths 
qi, ?2, qi- 

Its (squared) area can be expressed by 

qi Mqi + q2-qi) 



2 1 

\{qi + 92 - qi) 92 

= ^ + t;(o"i + cr 2 + cr ; )e + ^(cri<7 2 + a x ai + a 2 ai - ^)e 2 . (7) 

Expanding \J~&\ the series consists only of terms up to a 3 since af = 1. This suggests the 
following ansatz 

A t = c (e) + ci(e)(cri + o 2 + 07) + c 2 (e)(cr 1 cr 2 + o x Oi + o 2 Oi) + c 3 (e)cr 1 cr 2 (Ti , (8) 

where the coefficients c 4 - can be obtained by comparing with those in (7) (cf. Figure 2). 
Using (6) the measure (5) can be replaced by 

E exp[-m£ln(I + ^)], (9) 
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Fig. 2: Coefficients c 4 - obtained from first terms of series expansion (o) and semianalytically by 
comparison of coefficients (•). 



and the exponential in terms of e is reduced to 

J2H 1 + = ^i^o(e) 



m 



(10) 



with m 



mi 



m(e+ |e 3 ) 



+ 0(e 5 ) and iVi the total number of links. 



Inserting (8) and (9) into the partition function (4) yields 



Z 



Jexp{- J2(2Xc 1 + mi)ai - 
A ^[c 2 (cricr 2 + (J\(Ti + (J 2 (Ji) + c 3 cricr 2 cr ; ]} 



Jexp{- J2[(2\ci + mi)ai + 
+ ^\c 2 ((J 1 + cr 2 + cr 3 + (J 4 )(Ji + ^\c 3 ((J 1 (J 2 + cr 3 cr 4 )cr ; ]} 



11^ 



with J = exp( — AiV 2 c — Nim ) and iV 2 the total number of triangles. Thus, the path 
integral became the partition function of a system consisting of a spin 07 at each link /, with 
an external "magnetic field" and with two- and three-spin nearest neighbor interactions. 
Assigning the spin to the corresponding link of the original triangular lattice and drawing 
the interactions as lines a Kagome lattice is obtained (cf. Figure 3). Removing the term 
linear in a by a convenient choice of measure (mi = — 2\c\) and neglecting three-spin 
couplings we get the partition function of an Ising model on a Kagome lattice. To compute 
the critical coupling we abbreviate Q = — |c 2 A and transform the partition function Zk(Q) 
of the Kagome lattice to the partition function Zd(V) of the decorated honeycomb lattice 
by applying the star-triangle transformation [4]: 




Fig. 3: The triangular lattice is drawn by dashed lines and the Kagome lattice by solid lines. The 
spins are represented by o and • characterizing a typical configuration to be discussed in the text. 




Aexp[Q(cr 1 cr 2 + cr 2 cr 3 + <7 3 cri)] = ^ exp[Zcr(<7i + cr 2 + cr 3 )] . 

<T = ±1 
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For o"i = o 2 = cr 3 = 1 and <Ti = o 2 = — cr 3 = 1 we get 2cosh(3iv) = Aexp(3Q) and 
2coshZ/ = Aexp( — Q) } respectively. Any other choice of Oi coincides with one of these two 
cases. Hence, A and Q are determined from 



A 4 



+ 3) 5 



e 4Q = 2cosh(2iv) - 1 



(13) 



Carrying out this transformation at every vertex gives Zd(L) = A~ Nh Zk(Q) where Njj is 
the number of vertices of the honeycomb lattice. Furthermore, we can express Zd by the 
partition function Zh of the conventional honeycomb lattice via the decoration-iteration 
transformation [4]: 



L 



L 



K 



<7i a cr 2 (T\ cr 2 

^2 exp[Zcr(cri + cr 2 )] = /exp(AVi<7 2 ) . 



(14) 



(7 = ±1 



Similarly, we have 



I = 2e K 
> 2K = cosh(2iv) 



(15) 



and inserting (13) and (15) we obtain 



Z D (L) = [e^(e^ + 3p N »Z K (Q) = (2e K ) N »Z H (K) 
e 4Q = 2e 2K - 1 . 



(16) 
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Fig. 4: Spin expectation value as a function of the 2- and 3-spin couplings Q2 and Q3, respectively. 



The critical exponent of the honeycomb lattice is given by e 2Ac = 2 + y/3 and thus for the 
Kagome lattice 

e 4 ^ = 3 + 2^3 
Qc = -|(Ac 2 ) c = 0.4643 . (17) 

The critical coupling Q c is positive which reflects the fact that the antiferromagnetic 
Kagome lattice is always frustrated even at zero temperature. 

Switching on the 3-spin interactions we performed numerical simulations of the system 
in the coupling plane (Q 2 = — |Ac 2 ,(53 = ~~ 3^3)- Figure 4 shows the spin expectation 
value in a certain range of coupling constants. For Q 3 = the conventional Ising model 
on a Kagome lattice is recovered and for Q 3 ^ 0, Q2 < the 3-spin coupling removes the 
frustration leading to ordered phases: Depending on sign((5s) the spins denoted in Figure 

3 by open dots take the values a = ±1 whereas the full dots take a = =pl giving rise to an 
expectation value (cr) = ±|. 

We now return to the application of the spin system to the gravitational system. To 
investigate different types of measures the external field term has to be taken into account. 
Further the couplings depend on the value of e according to (8). Figure 5 depicts (cr) 
versus e for the uniform and the scale invariant measure and for the measure leading to a 
cancellation of the term linear in cr. The upper curve can be read off directly from Figure 

4 and represents only a small detail near the origin. The 3-spin coupling is in the entire 
range of e not strong enough to remove the frustration and therefore (cr) 0. Both other 
curves show a favored occurence of negative spins, corresponding via qi = 1 + eo\ to short 
links reflecting the expected tendency of the lattice to shrink [2]. 

The motivation for this exploratory study was to approximate quantum gravity by a 
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Fig. 5: (a) versus e for nil = —2c\ (no external field), m = (uniform measure) and m = 1 (scale 
invariant measure). In all cases the cosmological constant is A = 1. 



spin system. It is straightforward to extend this approach from the trivial two-dimensional 
Regge- Einstein action to scalar fields or physical spin fields. To allow for more than two link 
lengths an extension to Zjv is possible leading to an increase in the number of coefficients 
C{. Recently, in three dimensions the Regge-Ponzano model has been considered as a q- 
deformed su(2) spin network [5]. A generalization of our approach to higher dimensions is 
possible in principle but becomes more complicated. Notice that in two dimensions only 
the three links of a triangle are coupled in the action. In three dimensions one is faced with 
terms of 6 th order at least, and with additional contributions from the Regge action. The 
situation is even more complicated in four dimensions where one has to deal with 10 links 
in each simplex. 
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